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We propose an effect whereby an electric current along the interface between a ferromagnetic and 
normal metal leads to injection of pure spin current into the normal metal, if the magnetization- 
direction in the ferromagnet varies along the direction of current. For the specific example of a 
ferromagnetic domain wall, we compute the inverse spin-Hall effect voltage this spin current gives 
rise to when injected into a Pt layer. Furthermore, we show that this pure spin current leads to 
modification of the parameters that govern spin transfer and current-driven domain-wall motion, 
which can be use to optimize the latter in layered magnetic systems. This effect in principle enables 
control over the location of spin-current injection in devices. 

PACS numbers: 75.78.Fg 



Introduction — Spintronic devices make use of the 
spin degree of freedom to process and store information. 
Hence, the generation and detection of non-equilibrium 
spin accumulation and spin currents is of paramount im- 
portance. In particular, all-electric injection and control 
of spin currents at room temperature and without high 
magnetic fields is crucial for viable integration with and 
as extension of current technology To obtain spin 
currents, a large variety of physical mechanisms and ge- 
ometries are investigated jH [3] . One class of approaches 
relies on parametric pumping [3]. In these pumping ap- 
proaches a periodic (AC) excitation is transformed into a 
DC spin current. Examples are circularly polarized opti- 
cal photons [3l [SHZ] J magnons [SHIO] , acoustic waves [11] 
and single-domain ferromagnetic resonance [T2HT4] . 

In contrast to the pumping approaches described 
above, a spin accumulation can also be obtained via 
a static bias. A current through a ferromagnetic- 
nonmagnetic- metal (FM-NM) junction causes spin injec- 
tion into the nonmagnetic layer [15|, 116) . Room temper- 
ature injection in silicon was demonstrated in Ref. |17j . 
Spin-orbit coupling also opens the possibility to create 
spin currents using electric fields only. For example, the 
spin Hall effect generates a spin current transverse to a 
charge current [I8H2T] . 

The interplay between charge, spin and temperature, 
studied in the young field of spin caloritronics, yields 
novel ways to inject pure spin currents using temperature 
gradients. Thermal transport leads to the generation of 
spin accumulation via the spin Seebeck effect [12], spin 
dependent Seebeck effect [23] and spin Seebeck tunneling 

In this Letter we propose a mechanism for the injec- 
tion of a transverse spin current from a FM-NM interface 
by a charge current, as illustrated in Fig. [TJ Consider a 
(pinned) domain wall in the FM layer such that the direc- 
tion of magnetization m{x) depends on the coordinate 
along the wire, and such that an electric current flows 
along the wire. We denote the charge current density on 




FIG. 1. (Color online) A charge current along the interface 
between a ferromagnet containing a domain wall and a normal 
metal leads to the injection of a transverse pure spin current 
at the interface. The direction of the spin component of the 
injected spin current is in the direction of the helicity m x (j'c- 
V)m of the domain wall along the current direction, leading 
to a spin-accumulation in the normal metal. 



the interface by jc- It is the magnitude of this charge cur- 
rent density that governs the magnitude of the injected 
spin current. As we discuss in detail below, the change 
in orientation of the magnetization due to the presence 
of the domain wall will cause a transverse spin current 
ji^" into the NM given by 

f:'i^)-^^^i^)-^{'"s-'^)m{x), (1) 

where j™{x) is the injected spin-current density flow- 
ing perpendicular to the interface (and hence transverse 
to the charge current direction) with spin polarization 
in the direction of the helicity m x (jc • V)m of the 
magnetization texture along the current direction (in the 
y-direction for the situation in Figjl]). Here, mix) is 
the unit- vector magnetization direction of the FM-layer, 
Vg — —gLlJ-BPjc/'2Ms\e\ is the spin velocity, P is the spin 
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polarization in the ferromagnet, gL the Lande g-factor, 
/is is the Bohr magneton, is the saturation magneti- 
zation, and — |e| the charge of an electron. Furthermore, 
the parameter g — "qg^^ is governed by properties of the 
FM-NM interface, where g^^ is the real part of the mixing 
conductance in units fi~^m~^, and 77 = O (1) is a dimen- 
sionless parameter and Gq — 2e^ /h is the quantum of 
conductance. Eq. ([I]) is the main equation describing the 
effect we propose in this paper. 

In the next section derive an expression for the coeffi- 
cient g in terms of spin-dependent scattering properties 
of the FM-NM metal interface, which can in principle be 
evaluated for realistic interfaces. This derivation is some- 
what technical and some readers may wish to skip to the 
end of this section where we also give an estimate for 
g using a physical interpretation of the transverse spin 
current in terms of spin pumping. 

Scattering Theory — We consider a two-dimensional 
tight-binding model for the interface between ferromag- 
net and normal metal. The ferromagnetism is described 
with localized magnetic moments exchange coupled to 
the electrons, and connected on every site r = {ix, iy} to 
a metallic reservoir (i.e., the normal metal) with chemical 
potentials /ir. The system is described by the Hamilto- 
nian H = Hs + Hl + He representing electronic system 
in the FM, NM-leads and contacts between them. The 
first term is specified in terms of second-quantized oper- 
ators Cr^a{c\ a) that annihilate (create) an electron with 
spin a at site r 



ns = -Js 
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which describes nearest-neighbor (indicated by the brack- 
ets (•,•)) hopping with amplitude Js and coupling to the 
local magnetic moments rrip, where A is the exchange 
energy and t is the Pauli spin-matrix vector. The metal- 
lic contact is described by a set of one-dimensional leads 
at chemical potential /ir and modeled by the Hamilto- 
nian T-Ll = X^r^r where the Hamiltonian for lead-r is 



given by 



U'J")\o- 



where the 



and 



t 



are the 



hopping amplitude is Jl^ and ^ 

fermionic operators in the r"^ lead. Finally the contact 
between the ferromagnetic system and its leads is de- 
scribed by He = Z^r ■^ith 



at 3Lr 

-r,a"'aLr,<7 



(3) 



where dL^ denotes the last site of the lead and Jq is 
the amplitude of tunneling among both subsystems. An 
electric current flows through the FM by a difference in 
chemical potentials fiL,R that connect the left and right 
sides of the FM. 



The spin current flowing from the system to the leads 
is determined as the rate of change of spin density in the 
r-th site, i.e. ^ = 3^2 ^ E^,<t' (cLT-<T,<T'Cr,<T'), where a 
is the distance between sites. The general expression has 
been derived before 25j and is given by 



~dt 
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where N{e) — ^e'^/'^^T ^ ^ is the Fermi-Dirac distri- 
bution function with ksT the thermal energy, A'^(e) = 
G(+)(e)^r'"(e)G(-)(e) is the spectral-weight contribu- 
tion due to the lead at site r and the rate r'"(e) = 
j [l]'''^+)(e) — S''''-^^(e)] . Since we are assuming non- 
magnetic leads, the self-energy E'''*^+''(e) will be propor- 
tional to the identity in spin space. Its only non-zero 
matrix elements are HEl'l^^ ^, = -{Jc)'^e''''''^''^"'6a,a' /Jl,, 
with k'^{e)a = arccos [— e/2Ji^]. To carry out the explicit 
evaluation of Eq. Q in terms of the magnetization orien- 
tation rrij. it is convenient to decompose the Green's func- 
tions into spin-independent singlet and spin-dependent 
triplet parts, namely 



(5) 



Taking the trace over the spin indices in Eq. Q we dis- 
tinguish two contributions to the spin current, one com- 
ponent parallel to the magnetization vector and other 
transverse, denoted by js,r and respectively. The 
transverse spin current density induced by the magnetic 
texture to lead r is given by 



Js.r ^2 



de 



(27r) E i'^-f^" ) trr' (f) ("^r X 171^ 



(6) 

with the transmission probability ^^Jr'i^) = 

fiKA^)Gif'^\^)fiK'',r'{^)Gr'i~\^) for the spin- 
polarized part of the current flowing from the lead 
at site r to the lead at site r' through the FM. We 
now consider a zero net current (but nonzero spin 
current) flow into the leads, except for the left and right 
leads that have chemical potentials = £f + \e\V 
and = ep, respectively, with ep the Fermi en- 
ergy. At low temperatures and assuming the length 
scale of magnetization-orientation variation much 
greater than the inverse Fermi wavelength we see 
that in the continuum limit the only contributions 
to Eq. ^ are from neighboring leads for which 
rrir X TTir' m(x) X adxm{x). Keeping these contribu- 
tions, we find that the transverse spin-current density to 
lowest order in magnetization gradients satisfies Eq. ([l]) 
with g = a{l&Ms/gLl^BP) x (Got*(eF)/t(eF)), where 



(4) 
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t(ei?) is the total (i.e., summed for both spin channels) 
transmission probability. Both this transmission proba- 
bility, and the spin-dependent transmission probability 
are taken between leads at neighboring site and are 
taken at the Fermi energy in the low-temperature limit. 
Also note that these transmission probabilities are 
evaluated for the homogeneous ferromagnetic state, as 
Eq.(|6]) is already first order in magnetization gradient. 

We obtain an estimate for the coefficient g by con- 
necting the transverse spin current [c.f. Eq. ([T])] to spin 
pumping by a precessing magnetization. To make this 
connection we consider a helical magnetization mh_e\{qx+ 
Lot) precessing with frequency u. We assume the wave- 
length of the helix to be much larger than the spin dif- 
fusion length in the ferromagnet, qXgd ^ Ij such that we 
can use the spin-pumping expression for a single-domain 
ferromagnet [T21 US] 



hg^^ dm 



(7) 



for the spin current pumped through the FM-NM 

interface by a time-dependent magnetization. Now we 
make a Galilean transformation to the frame moving with 
velocity uj/q. In this frame we have a static magnetiza- 
tion and a spin velocity of —lu / q and obtain Eq. ([!]) with 
g/g^^ = 1 by demanding that the injected spin current 
does not change. Note that Galilean invariancc is in gen- 
eral broken due to disorder leading to a g/g^^ 7^ 1 but 
of order unity. From this argument we can also make an 
estimate for the magnitude of the transverse spin cur- 
rent injected by a domain wall of width Adw For narrow 
domain walls with Adw of the order of a few nm and cur- 
rents for which Vs is of the order of 100 ni/s (that can be 
obtained in experiments), the injected spin current den- 
sity would be equivalent to pumping with a ferromagnetic 
resonance frequency of Ws/Adw ~ 100 GHz. 

Inverse Spin Hall Effect — We now consider how 
the transverse spin current can be probed by the in- 
verse spin Hall efi'ect (ISHE). We consider a pinned 
domain wall with a width Adw in a perpendicular 
magnetized FM, as shown in Fig. [T] where iy9dw 
is the angle of the magnetic moments in the wall 
with the X-axis. We model the magnetization by 
tTTidwix) = (cos iy9dw sin0, sin (/?dw sin0, cos6')"^ , where 6* = 

2arctanexp ^Q^xf^)' ^"^^ Q charge of the do- 

main wall. The injected spin-current will result in a 
nonzero spin accumulation fisix,z) which can be found 
by solving the spin-diffusion equation in the NM [13], 
which in the static limit r^fdiJis/dt <^ 1, is given by 



\2 ' 



(8) 



where Agd = Vl^s^sf is the spin-diffusion length in the 
NM, and and Tgf are its spin-diffusion constant and 
spin-flip time, respectively. The boundary conditions for 



Eq. ([8]) enforce continuity for the spin current and are 
given by 



dzfJ'Ax,z)\^^^^ = 0; 
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(9) 
(10) 
(11) 



where Lat, (In and a are the length, thickness and con- 
ductivity of the NM respectively, js°zi^) is the net spin- 
current into the NM. The solution for the spin accumu- 
lation yields fis{x,z) — J dx'K{x — x' , z)j^°^{x') where 
the fourier transform of the kernel K{x, z) is given by 



k{k^,z) = - 



GqX, 



sd 
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(12) 

The total spin current is the sum of the spin-injection and 
spin-pumping contributions, as given in Eqs. ([1]) and 0, 
and a spin-current j^'^'^^ m the opposite direction due to 
the induced spin- accumulation on the NM side of the 
interface. Thus = jl 



.-pump 



^•back 

Js,z ' 



where the 



back-flow of spins trough the interface is given bv[26] 



•back 
J s^z 
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Hsix,z = 0), 



(13) 



where we neglected the imaginary part of the mixing 
conductance. This equation is justified for realistic 
interfaces [22]. Using Eq. (12) in Eq. (13) we get an 



expression for the total spin-current in terms of the spin- 
injection and pumping contributions only. Note that here 
the pumping contribution is zero as we consider a static 
pinned domain wall. We come back to this contribu- 
tion when we consider the interfacial enhancement of spin 
transfer torques. 

The ISHE [351 HH] gives a voltage signal from a spin 
current via jl^^^ = x ^SHe'''"i"j, with the spin current 
given by jfj — ^dfi"/dx^ and 6'sh the spin-Hall angle. 
We average the voltage generated via the ISHE over the 
thickness djsf of the NM-layer. The voltage difference due 
to the ISHE in the x-direction is then 



AyjSHE 



2d 



N\e\ 



T^SH 



(cosh ^ 



l)Qcos (/3dv 



sinh ^ 



cosh ■ 



which holds for the domain wall far from the edges of 
the NM-layer. Note the thickness of the Pt-layer in- 
fluences the signal considerably [30 , but that the result 
does not depend on the width of the domain wall be- 
cause the ISHE induced voltage is in the x-direction. 
For a Co-Pt interface we take gl'^ « 4 • 10^5r?-^m-^[3T]. 
CTpt = 9.5 • l{f'Q.^^in~^ , Asd = 1.5nm, Adw ~ 5nm and a 
platinum thickness of d^ = 3nm, ffgu = .05 and obtain 
AVjf^^^/Ty ~ 40nV. For a current density in the cobalt 
of « 10^2 Am"^ 
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Interfacial enhancement of spin transfer — Another 
way to observe the injection of the spin current is the 
effect it has on the motion of the domain wall. This is 
important for the understanding of domain-wall dynam- 
ics in layered magnetic materials, that are the subject 
of ongoing research [32] ■ The total spin current ejected 
from the FM-layer reduces the angular momentum of the 
FM, inducing a torque which modifies the domain-wall 
dynamics. The Landau-Lifshitz-Gilbert equation in the 
presence of this additional torque is given by 



dm 

~dt 



dm 

dx 



ao^M X 



l3oVsm X 



dm 
dm 



^m X ^^^^ 
h dm 



1 



dx IvLd, 



•net 

Js,z 7 



(14) 



where ao and /3o are the bulk Gilbert damping and non- 
adiabaticity parameter respectively. The FM has thick- 
ness dp and hard-axis anisotropy K± which we take to 
be along the y-axis and is included in the energy func- 
tional -Emm that also contains the exchange and easy 
axis anisotropy that set the domain-wall width. From 



the LLG equation [Eq. (14)] we obtain the equations of 
motion for the collective coordinates of the domain wall, 
which are given by 



dw 



~2h 



sin 2(^dv 



Adw 



(15) 
(16) 



where a^.i^ and are given by 



ao 



dxdx n(a; — x')j 



'\ .-in+pump 



{x) 



-{x ); 



5^ 



(17) 



with ^ = {r, Lp\ and where n(a; — x') is given by 



I{{x - x') 



27r 



47rGo 



Jk^(x~x') 



Note that the integrals T^^r are functions of the dimen- 
sionless parameters Adw/Asd, dj^/Kd and 5^"^Asd/47rcr. 
The average velocity of the current-driven domain wall 
is given by 



sign 



1 



13^ 



Re 



/3. 



where the critical velocity is given by Vc — K±Xdw/'2fi- 
The ratio determines the qualitative behavior of 

the domain- wall velocity as a function of current. In Fig. 
[2] this ratio is shown as a function of the thickness of 
the normal metal layers in the multilayer. Note that in 
the limit of a large ratio Adw/Asd our result for alpha co- 
incides with interfacial enhancement of Gilbert damping 
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FIG. 2. (Color online) The ratio P^/a^ is shown as a function 
of the thickness c^n of the normal metal layers in units of its 
spin diffusion length Asd ~ 250nm, for a Cu-Py-Cu multilayer. 
The head-to-head domain wall in the Py has a width of 100 
nm. For g/g^^ = 2, where g^^ = 1.6 ■ W^'^Q'^m'^ is the 
mixing conductance of a Py-Cu interface the ratio, the ratio 
Ptp/ctif = 2 coincides with the bulk ratio with ao = 0.006. 



in single-domain magnets obtained by Tserkovnyak et al. 

[una. 

Discussion and Conclusions — In this Letter we de- 
scribed a mechanism for transverse spin injection across 
a FM-NM interface that is induced by a magnetization 
gradient along a current, and focused on the example of 
a pinned domain wall. We note here that if the domain 
wall is not pinned and moves with velocity Wdw, the in- 
jected spin current is altered which is described by the 
replacement Vs ^ Vs — t'dw hi Eq. ([T]). The position 
of the injection is dependent on the position of the do- 
main wall which could lead to controllable local spin in- 
jection. Moreover, due to the spin accumulation induced 
in the NM layer (s) the measured a and /3 parameters for 
domain- wall motion for multilayer systems will be differ- 
ent than bulk values. We expect that this effect plays a 
role in all thin- film measurements of domain- wall dynam- 
ics. In fact, large values of beta are typically reported in 
such systems [33] , pointing to the possibility of interfacial 
enhancement. 

Other spin-injection mechanisms, like the spin- 
dependent Seebeck effect [23] or diffusive spin 
injection [16], typically induce a spin current in the 
NM-layer with the spin direction parallel to the magne- 
tization in the ferromagnet. Therefore the spin-injection 
we discussed in this Letter, which induces a spin 
current polarized in the direction of the helicity, is 
distinguishable. 

We also note that in principle there is a spin current 
with spin polarization in the direction of {vs • V)m that 
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we have ignored as it oscillates and averages out to small 
values when integrated over position, and is determined 
approximately by the imaginary part of the mixing con- 
ductance which is small for realistic interfaces. 

The domain-wall induced spin injection could be used 
in a memory device consisting of a FMi— NM — FM2 tri- 
layer where FMi is a ferromagnet with a domain wall 
adjacent to FM2, which is a monodomain. For the geom- 
etry as shown in Fig. [T]the injected spin-current has the 
same symmetry as the SHE. We can estimate the magni- 
tude of the effective spin Hall angle for the spin-injection 
to be 6* = |2ejs//ijc| ~ 0.008 for a Co-Pt interface, where 
the spin current was estimated by replacing the gradient 
in Eq. ([T]) by 1/ Adw • The presence of this spin current 
in addition to the ISHE spin current could be used to 
control the switching of the single domain FM2 by the 
presence of a domain wall in FMi. 

The reverse process of spin-injection occurs as well. A 
spin current flowing with a transverse spin direction into 
a ferromagnet at the position of the domain wall will lead 
to a voltage difference over the ferromagnetic strip. In 
this way domain walls could be used as a local moveable 
sensors of spin current. In future work we plan to inves- 
tigate transverse pure spin currents associated with heat 
currentsin the same geometry as discussed here. 
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